Singularities and the distribution of density in the Burgers/adhesion model 
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We are interested in the tail behavior of the pdf of mass density within the one and d-dimensional 
Burgers/adhesion model used, e.g., to model the formation of large-scale structures in the Universe 
after baryon-photon decoupling. We show that large densities are localized near "kurtoparabolic" 
singularities residing on space-time manifolds of codimension two (d < 2) or higher (d > 3). For 
smooth initial conditions, such singularities are obtained from the convex hull of the Lagrangian po- 
tential (the initial velocity potential minus a parabolic term) . The singularities contribute universal 
power-law tails to the density pdf when the initial conditions are random. In one dimension the 
singularities are preshocks (nascent shocks), whereas in two and three dimensions they persist in 
time and correspond to boundaries of shocks; in all cases the corresponding density pdf has the 
exponent —7/2, originally proposed by E, Khanin, Mazel and Sinai (1997 Phys. Rev. Lett. 78, 1904) 
for the pdf of velocity gradients in one-dimensional forced Burgers turbulence. We also briefly con- 
sider models permitting particle crossings and thus multi-stream solutions, such as the Zel'dovich 
approximation and the (Jeans)-Vlasov-Poisson equation with single-stream initial data: they have 
singularities of codimension one, yielding power-law tails with exponent —3. 
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sopra quei gruppi, sopra quegli ammassi, 
quel seminio, quel balenio di stelle 

Giovanni Pascoli, from La Vertigine 



I. INTRODUCTION 

In 1970 Zel'dovich introduced a simple model for ex- 
plaining features of the nonlinear formation of large-scale 
structures in the Universe Just after the baryon- 
photon decoupling in the early Universe, there may have 
been a rarefied medium formed by collisionless dustlikc 
particles without pressure, interacting only via Newto- 
nian gravity Q . The appropriate mathematical descrip- 
tion, the equation for a self-gravitating gas in an expand- 
ing three-dimensional universe, has so far been studied 
mostly by numerical simulations. The Zel'dovich ap- 
proximation, to which we shall return in Section M, is 
far simpler and involves basically particles moving in 
straight lines, just as rays in geometrical optics (see also 
Refs. ^,|^1). As a consequence caustics are formed, the 
simplest of which are pancakes, near which the mass den- 
sity is very large. Zel'dovich, Arnold and their collabo- 
rators were mostly interested in the nature of the singu- 
larities resulting from the model and classified them us- 
ing catastrophe theory and Lagrangian singularity the- 
ory ||||-^|. Kofman et al. ^ studied the probabilistic 
aspects and determined the probability density function 
(pdf) p(p) of the density of matter. For Gaussian initial 
fluctuations (at decoupling) they obtained a p~ 3 law for 
the tail at large densities. They also gave a simple heuris- 
tic argument relating the law to the divergence with 



exponent —1/2 near the pancakes. 

One difficulty with the Zel'dovich approximation is 
that the pancake structures, which are formed after the 
first particle crossing occurs, rapidly smear out (see Sec- 
tion IV), whereas in reality massive pancake-like struc- 
tures are found to be quite long-lived. The gravita- 
tional dynamics of pancakes is indeed incorrectly cap- 
tured within the Zel'dovich approximation (see, e.g., 
Rcf. jl0|]). This has led to the introduction of the ad- 
hesion model of Gurbatov and Saichev pd| , p^ ] in which 
particle, upon crossing, stick together (adhere). The ad- 
hesion model is just the multi-dimensional Burgers equa- 
tion, taken in the inviscid limit v — > 



d t v + (v V) v = z/V 2 v 
v = -VV>, 



d tP ^ + V 



0. 



(1) 
(2) 
(3) 



Here, v = (u, v, w) is the velocity, if) the (velocity) po- 
tential and j(/ E ) the Eulerian (mass) density (the initial 
Lagrangian density, which is quasi-uniform in the cos- 
mological problem, is denoted po). As is well known, 
the Burgers equation in the limit v — > produces shocks 
along surfaces (in three dimensions) on which the density 
is infinite and across which the velocity is discontinuous. 

The question we intend to address is the behavior at 
large p's of the pdf of the density when using the one- 
and multi-dimensional Burgers/adhesion model with ran- 
dom and smooth initial conditions, having, e.g., Gaus- 
sian statistics with an exponentially decreasing spectrum. 
(The case of non-smooth (e.g. Brownian) initial condi- 
tions has been considered in Refs. |l^-[n|.) In one di- 
mensional decaying Burgers turbulence, the density and 
the Eulerian velocity gradient have a simple relation (cf. 
(|l6|)). Hence, the density pdf is deducible from the pdf 
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of the velocity gradient. It was shown in Ref. |16|] that 
the latter has a tail at large negative values which is 
a power law with exponent —7/2. Actually, this law 
has been first proposed in Ref. |17j] for randomly forced 
Burgers turbulence. For the forced case, the functional 
form of the pdf has been the subject of considerable 
controversy and the question is not yet completely re- 
solved (see Rcfs. |l7] ^8|). In particular E and Vanden 
Eijnden [^l],^2| developed a probabilistic formalism that 
copes with the delicate problems arising in the limit of 
vanishing viscosity when shocks are present, proved that 
a < — 3 and made a good case for a = —7/2. 

The paper is organized as follows. In Section || we 
present general background material about the multi- 
dimensional Burgers equation and show why, paradox- 
ically, shocks are not generally responsible for large (but 
finite) densities, which generically arise from other types 
of singularities. In Sections [ill] and IV we present the 
one-dimensional and the multi-dimensional case, respec- 
tively. In the final Section ^ we compare the predictions 
of various models used in the cosmological literature: the 
Burgers/adhesion model, the Zel'dovich model and the 
(Jeans)-Vlasov-Poisson model. 



II. THE SOLUTION OF THE BURGERS 
EQUATION AND ITS GEOMETRICAL 
CONSTRUCTION 

We shall here be exclusively interested in the solu- 
tion to the d-dimensional Burgers equation (|l|)-(^), in the 
limit v — > 0, with given initial data vo(x) = v(x, 0) = 
-V^o(x) and po = p( E )(x,0). 

We begin with the velocity. The Hopf |29[ and Cole 
p0| transformation allows an explicit integral represen- 
tation of the solution for v > 0. By steepest descent, a 
well known "maximum representation" is obtained in the 
limit v ->■ iHJl 



V>(x, t) = max 
q 



^o(q) 



(x 



\2l 



2t 



(4) 



It is easily seen that the point q at which the maximum 
is achieved is the Lagrangian point associated to the (Eu- 
lerian) point x at time t. Indeed, by setting the gradient 
of the r.h.s. of (0) to zero and using (0), we obtain 



x = q + fv (q). 



(5) 



In other words, x is the position at time t of the fluid 
particle starting at q and retaining its initial velocity 
vo(q). 

The problem is that (^) is valid only for regular La- 
grangian points, that is points which have not been cap- 
tured by a shock by time t. There is another construction 
of the solution which brings out the geometrical nature 



of the problem with shocks. Let us define the Lagrangian 
potential 



v(q,f) = — +^o(q) 



It follows from (W) that 



— = max [<p(q, t) 
2 q 



•x-qj 



(6) 



(7) 



The r.h.s. of (|7|) is seen to be the Legendre transform of 
the Lagrangian potential Q . 

An important property of the Legendre transformation 
is that the r.h.s. of is unaffected if we replace the 
Lagrangian potential <^(q, t) by its convex hull <£> c (q, i). 
The Eulerian singularities of the solution are determined 
by the structure of the convex hull. 

Since we shall make extensive use of convex hulls in 
this paper, let us define the matter precisely. Let g(q) 
be a real function of q defined over a convex domain T> 
of IR'' (for example, the whole space). We say that <?(q) 
is convex if, for any q G T>, q' G V and < 8 < 1, we 
have 



g (0q! + (1 - 0)q) > 0g(q[) + (1 - %(q). 



(8) 



Let y(q) be an arbitrary function. We define its convex 
hull <£> c (q) as 



Vc(q) = min5(q) 



(9) 



the minimum being taken over all functions g(-) > ip{-) 
which are convex. 

In other words, the graph of <p c (q) is obtained by 
tightly pulling a string (in one dimension) or an elas- 
tic sheet (in two dimensions) around the graph of y(q) 
(see Figs, [l] and ||). 
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FIG. 1. Lagrangian potential and its convex hull in one 
dimension. The graph of the convex hull contains regular 
parts of the Lagrangian potential and segments touching the 
original graph at two points, lying over shock intervals. 
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FIG. 2. The convex-hull construction in two dimensions 
(adapted from Ref. [EHl); the graph of the convex hull con- 
tains regular parts of the Lagrangian potential (corresponding 
to regular points in the Eulerian space), pieces of ruled sur- 
faces (corresponding to shock lines), triangles (corresponding 
to shock nodes) and "kurtoparabolic points" (corresponding 
to termination points of shock lines). 

In one dimension the graph of the convex hull con- 
tains (i) parts of the graph of <p(q) and (ii) segments 
touching the original graph at two points. The former 
correspond to regular (Lagrangian) points and the latter 
to points having fallen into a shock. In two dimensions 
the convex hull consists generically of four kind of ob- 
jects: (i) parts of the original graph, (ii) pieces of ruled 
surfaces, (iii) "kurtoparabolic points" , to which we shall 
come back, and (iv) triangles (see Fig. ||). The associ- 
ated Eulerian objects are, respectively, (i) regular points, 
(ii) shock lines, (iii) end points of shocks and (iv) shock 
nodes. For the three-dimensional case, see Ref. |34j and 
references therein. 

A more complete description of singularities is ob- 
tained by considering the metamorphoses of singularities 
as time elapses. A complete classification in two and 
three dimensions may be found in the appendix (supple- 
ment 2) by V.I. Arnold, Yu.M. Baryshnikov and I. A. Bo- 
gayevski of Ref. |L2| . 

Let us just observe at this stage that, when the initial 
potential ipo(<l) is a sufficiently smooth function of q (as 
we assume in this paper) , the Lagrangian potential given 
by (^) is necessarily convex when t is sufficiently small. 
Indeed, the Hessian (determinant of the Hessian matrix 
of second space derivatives) 



H(q, t) ee det 



aV(g,*) 

dqidqj 



t d det 



gVo(q) 
dqidqj 

(10) 



remains very close to its initial value (— l) d for short 
times. As long as the Hessian does not change sign, 
the convexity is the same as for the initial paraboloid 



^(q,0) = -|q| 2 /2. 

We also introduce two Lagrangian maps from q to x. 
The naive Lagrangian map L t is just given by (||). The 
(proper) Lagrangian map is 



L t : q h-> -V q <^c(q, t) = x(q, t). 



(11) 



At regular points, where ip(q,t) and ip c (q,t) coincide, so 
do the two Lagrangian maps. If, however, q is not a reg- 
ular point, the Lagrangian map transforms it into the ap- 
propriate Eulerian shock location, whose determination 
requires the knowledge of the convex hull, that is, a global 
geometrical construction. The map q i— > x = £ t q is in- 
vertible only at Eulerian points which are not on shocks 
(otherwise there is more than one Lagrangian point which 
is mapped into x) . The Jacobian of the Lagrangian map 
at regular points is defined as 



J(q, t) ee det 



dxi 



(12) 



It follows from ( |lfj| ) and ( |TT| ) that it is just (— l) d times 
the Hessian -ff (q, t) of the Lagrangian potential. 

We turn now to the determination of the (Eulerian) 
density p( E )(x, t). Mass conservation implies that, if x is 
not on a shock, 



(X,f) 



Po 



J(q,i) 



(13) 



If x is on a shock, the density is of course infinite. This 
does not, however, imply that large but finite densities 
are generally obtained near shocks. Indeed, for p^ E ^(x, t) 
to be large, the Jacobian J(£ t " 1 x, t) and thus the Hes- 
sian H(C^ic, t) must be small. At any given time t, 
this happens only near the (d — l)-dimensional mani- 
fold of vanishing Hessian. Arbitrarily close to such a 
"parabolic" point there are generically hyperbolic points 
where the surface defined by <p(q) crosses its tangent (hy- 
per)plane and which, therefore, do not belong to its con- 
vex hull. As we shall see in Sections HI and |y| there 
can be "kurtoparabolic points" with vanishing Hessian 
which are at the boundary of regular regions. (In Greek, 
Kvproa means "convex" ; hence the proposed name. In 
Ref. |}4| they are called ^-points.) In one dimension 
kurtoparabolic points correspond to preshocks |35|] and 
exist only at discrete times; in two and three dimensions 
they generally persist for a finite time and are associated, 
in the Eulerian space, to boundaries of shocks (termina- 
tion points in two dimensions and edges in three dimen- 
sions). Large densities are obtained exclusively in the 
neighborhood of kurtoparabolic points. 



III. THE ONE-DIMENSIONAL CASE 
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A. Preshocks 

In one dimension, following fluid dynamical tradition, 
we denote the Lagrangian coordinate by a and the ve- 
locity by u. We denote by uo(a) — —dipo(a) / da the ini- 
tial velocity, which is assumed to be random and suffi- 
ciently smooth. It is convenient, but not essential, to 
assume periodic boundary conditions, homogeneity and 
a vanishing mean velocity. We denote by po the initial 
background density, taken deterministic and uniform. At 
regular points (outside of shocks) the Eulerian velocity 
and density are given implicitly by 

u(x,t) = u (a), p W( X)t ) = J^, (14) 
x = a + tu (a). (15) 

Since d x u = (d a uo) (d a x) 1 , it follows that the Eulerian 
density can be written in terms of the Eulerian velocity 
gradient d x u O] 



p^(x,t) =p a (l-td x u(x,t)). 



(16) 



In Ref. (T(| it was shown that the pdf of the Eulerian 
velocity gradient has a -7/2 power law at large negative 
values. Hence, the pdf p(p) of p^ has also a -7/2 power 
law, but at large positive values. The proof given in 
Ref. |l6| was rather detailed. Here, we give a simplified 
derivation, adapted to the case of the density. Further- 
more, we shall work mostly with the potential and nor- 
mal forms near singularities, to prepare the ground for 
the multi-dimensional case. 

From @ and |t|) we have, at regular points, 



pW(L t a,t) = 



Po 



l-td 2 il>o(a)/da 2 



(17) 



Large values of p^ are thus obtained in the neigh- 
borhood of Lagrangian points with vanishing Jacobian, 
where d 2 ipo(a)/da 2 = 1/t. Once mature shocks have 
formed, the Lagrangian points with vanishing Jacobian 
are inside shock intervals and thus not regular. The 
only kurtoparabolic points (points with vanishing Jaco- 
bian at the boundary of regular regions) are obtained at 
preshocks, that is when a new shock is just born at some 
time t*. Preshocks, play a central role in the -7/2 law of 
E, Khanin, Mazel and Sinai for the forced case jL7). Such 
points, denoted by a*, are local negative minima of the 
initial velocity gradient, characterized by the following 
relations 



da 2 ' 



d 3 iP 
da 3 



(a*) - 0, 



dSh 
da 4 



(a*) < 0. (18) 



There is however an additional global regularity condi- 
tion that the preshock point a» has not been captured 
before i» by a mature shock. This may be written in 



terms of the convex hull ip c of the Lagrangian potential 

ip, as 



(p(a*,t») = —al/2 + #0(0*) = ip c {a*,U). 



(19) 



As shown in [[l6|, this global condition affects only con- 
stants and not the scaling properties of p{p) at large p's. 

We can now Taylor expand the Lagrangian potential 
and the Lagrangian map near the space-time location 
(a*, t*). By adding a suitable constant to the initial po- 
tential, by performing a suitable translation and also a 
Galilean transformation canceling the initial velocity at 
a*, we may assume that a* = ^>o(a*) = difto(a*) j da = 0. 
We then obtain, to the relevant leading order, the follow- 
ing "preshock normal forms" 



<p(a,t) 

x(a, t) 
J(a,t) 



2 

-ra 

-T - 



-4Ca 3 , 



where 



(20) 

(21) 
(22) 



(23) 



The Lagrangian potential, together with its convex hull, 
are shown in Fig. |[ 

<P(a,t) 
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FIG. 3. Normal form of the Lagrangian potential in the 
neighborhood of a preshock in one dimension. At the time 
of the preshock (r = [t — t*)/t* = 0), the Lagrangian po- 
tential changes from a single extremum to three extrema and 
develops a non-trivial convex hull (shown as a dashed line). 

Note that at t = i* there is a degenerate maximum 
with quartic behavior and that, immediately after t*, for 
r > 0, convexity is lost and a shock interval is born. 
Given the symmetry, resulting from our choice of coordi- 
nates, the convex hull contains a horizontal segment ex- 
tending between the two maxima a± = ±(— r/(4^)) 1 / 2 . 
Note that for r > the Jacobian vanishes at two lo- 
cations ±(— t/(12(')) 1 / 2 which are within the shock in- 
terval and are therefore irrelevant as far as the Burg- 
ers/adhesion model is concerned (although they become 
relevant when particle crossing is permitted; see Sec- 
tion 0). 

From ( p2| ) we see that the density po/J has a a~ 2 singu- 
larity in Lagrangian coordinates at t = t* (r = 0). Since, 
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by (Elf), the relation between a and x is cubic at r = 0, 



the density p^\x,t*) 



oc |x| -2 / 3 which is unbounded. 
For any t 7^ i* the density remains bounded, except 
at the shock location. For t < 0, this follows immedi- 
ately from (22), which implies p^ < po/\ T \- For t > 0, 
the exclusion of the shock interval requires \a\ > a+. 
Hence, p^ < po/(2r). It is clear that large densities 
are obtained only in the immediate neighborhood of the 
preshock. More precisely, it follows from (21) and (|2T 
that p^ > p requires simultaneously 



|r| < ^ and \x\ < (-12C)~ 1/2 ^ 
P \P 



3/2 



(24) 



which become very small intervals around the spatio- 
temporal location of preshocks when p is large. 



B. The -7/2 law in one dimension 

So far, we have looked at the question of large densities 
from a deterministic point of view. We turn to the prob- 
ability to have a large density at a given Eulerian point 
x and a given time t. We shall calculate the cumulative 
distribution 



P>(p;x,t) = Prob {p( E )(x,t) > pj 



(25) 



from which we obtain the pdf p(p; x, t) — —d p P > {p\ x, t). 
In the random case each preshock has a random Eulerian 
location , occurs at a random time and has a random 
C < coefficient (there is also a random velocity u* of 
the preshock but this is easily seen to be irrelevant for 
our purposes). Only those realizations such that x* and 
are sufficiently close to x and t will contribute large 
densities. Denoting by ps(x*,t*, £) the joint pdf of the 
three arguments, which is understood to vanish unless 
£ < 0, we have 



P>(p;x,t) 



p( E '>(x,t)>p 



P3(x*,t*, C) dx* dt* dC- (26) 



(If homogeneity is assumed pz does not depend on x„; 
the case of homogeneity extending over the whole space 
can be obtained by letting the spatial period L — > oo; 
we must then also replace p?, by nlt^Q/L where n is a 
number of preshocks per unit length rather than a proba- 
bility; similarly, P > is then a probability per unit length.) 
Because of the very sharp localization near preshocks im- 
plied by (l24|), for large p's, we may replace ps (x* , i» , £) 
by ps{x, t, £). Using then, in a suitable frame, the normal 
forms (|2^)-(|2^) we can rewrite (|26| ) as an integral over 
local Lagrangian variables a and r and obtain 



P > (p;x,t)^ / t(-r- 12(a z )p 3 {x, t,QdadTd(. (27) 
Jd 

Here, the domain D is the set of (a, t, C) such that 



-4C 



< a" < 



-12C 



Po 
P 



(28) 



The right part of (28) expresses that the density exceeds 
the value p, while the left part (which is trivial when 
r < since £ < 0) excludes the shock interval ]a_,a+[. 
In (27) the factor — r — 12£a 2 is a Jacobian stemming 
from the change to Lagrangian space variables and the 
factor t stems from the change of temporal variables. The 
integration over a and r can be carried out explicitly, 
yielding 



P>(p;x,t)~C(x,t) 



5/2 



(29) 



C(x,t)=At \(\- 1 ^ P3 (x,t 1 C)dC (30) 



where A is a positive numerical constant. Thus, for any 
x and t, the cumulative probability of the density follows 
a p -5 / 2 law. Hence, p(p;x,t) oc p~ 7 ^ 2 , as p — > 00, which 
establishes the —7/2 law for the pdf. Note that, contrary 
to the derivation in Ref. |l6f| , we did not use homogene- 
ity. With this additional assumption, p% and thus C{x, t) 
become independent of x. 

Taking into account higher-order singularities does not 
influence this result. Consier, e.g., quintic-root preshocks 
arising from degenerate inflection points in the initial ve- 
locity, at which the second, third and fourth space deriva- 
tives all vanish. In one dimension such singularities are 
not generic in the deterministic case but could neverthe- 
less contribute in the random happens in Berry's 
"battle of catastrophes" (3(| . The exact vanishing of two 
more derivatives has probability zero but there is a finite 
probability that the third and fourth velocity derivatives 
have values small enough to give the preshock an approxi- 
mately quintic-root structure. We found that such events 
contribute only subdominant corrections to the -7/2 law. 



IV. TWO DIMENSIONS AND BEYOND 

A. Preshocks in two dimensions 

In two dimensions we use the notation q = (a, b) and 
x = (x, y) for Lagrangian and Eulerian coordinates. It 
follows from (|To| ) that the first singularity happens at 
the time which is the inverse of the largest positive 
eigenvalue of the Hessian matrix of the initial potential 
rpo(a,b). No generality is lost by making the following 
assumptions: (i) the maximum is achieved at the origin, 
(ii) the potential and its gradient vanish at the origin, (hi) 
the maximum eigenvalue is equal to one (i, — 1) and (iv) 
the eigendirections of the matrix of second derivatives are 
the a-axis for the eigenvalue 1 and the 6-axis for the other 
eigenvalue 1 — \i with p > 0. Using this, we can Taylor 
expand the initial potential to the relevant (fourth) order: 
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-0o(a, b) ~ — + (1 - /x) y + ™ 3 + /?a 2 fr + 7^ 2 + Sb 3 
+Ca 4 + 7?a 3 & + 6a 2 b 2 + nab 3 + pb 4 . (31) 

Expressing that the matrix of second derivatives has its 
largest eigenvalue at a — b — 0, we find that a — (3 — 
and that the quadratic form 

Q(a, b) = 12^Ca 2 + + (2p9 + 4 7 2 )6 2 (32) 

must be definite negative, thereby putting certain restric- 
tions on p, C, 9 and 7 which will henceforth be assumed. 

We can now write the corresponding normal form of 
the Lagrangian potential tipo(a,b) — (a 2 + b 2 )/2, at time 
t = 1+t for small t, which includes all the relevant terms 

a 2 b 2 

(p(a, b, t) ~ Ty - p— + 7a6 2 + Ca 4 + r/a 3 b + 6a 2 b 2 . 

(33) 

For r < the surface defined by the Lagrangian potential 
has a single maximum at a = 6 = 0. At r = this 
maximum still exists but is quartically degenerated in 
the a-direction. For r > 0, the origin turns into a saddle 
and two new maxima appear (for values of a = O (t 1 ^ 2 ) 
and of b = O (t 3 / 2 )). The general aspect of the surface 
is shown in Fig. It is clearly not convex; hence, a new 
shock is born. 




FIG. 4. The Lagrangian potential just after a preshock in 
two dimensions. The continuous line is the separatrix between 
the regular part and the ruled surface of the convex hull. The 
dotted line corresponds to the vanishing of Jacobian of the 
Lagrangian map. They touch at two "kurtoparabolic" points 
A and A'. 

The Hessian, now equal to the Jacobian of the La- 
grangian map, which controls the density at regular 
points is given by 

H(a, b, t) ~ —rp — Q(a, 6), (34) 

where Q is given by (|32|). 



Consider first the situation at t = t*. The Hessian van- 
ishes and, hence, the density becomes infinite at the ori- 
gin. This is the singularity called A 3 by Arnold P,p|,^3[ , 
for which the mean density in a small disk of radius r 
around the origin (in the the Eulerian space) is easily 
shown to diverge as p(r) oc r -2 / 3 . 

Shortly after i*, if we (incorrectly) use y(a, b, t) rather 
than its convex hull, we find that the density becomes 
infinite in Lagrangian coordinates along a zero-Hessian 
curve of approximately elliptical shape. On the surface 
defined by the Lagrangian potential ip(a,b,t), the corre- 
sponding points are parabolic (they are shown as a dot- 
ted line on Fig. |]). When approaching such a line of 
parabolic points, one has an Ai singularity in the sense of 
Arnold, for which p(r) oc r -1 / 2 . Actually, nearly all the 
points on this line are "hidden under the convex hull" . 
Constructing the convex hull of ip(a, b, t) is not a local 
operation and thus, in general, not elementary. Let us 
just illustrate what can happen in the simpler case where 
7 = 77 = 0, which has an additional symmetry. We then 
have 

<p(a,b,t)c£T^-Hj +(a 4 + ea 2 b 2 , (35) 

which is even in both a and b. The conditions of nega- 
tive dcfinitcncss of the quadratic form Q arc then ( < 
and 9 < 0. It follows from the symmetries and the fact 
that lines of constant a are parabolas that the convex 
hull contains a piece of ruled surface made of segments 
parallel to the a-axis. These segments are connecting the 
two maxima of sections at constant 6, which exist for any 
b 2 < —t/(29). The horizontal projections of these end 
points define the separatrix between regular Lagrangian 
points and points absorbed into the newly created shock. 
It is the ellipse 

r + 4Ca 2 + 26>6 2 = 0, (36) 

obtained by requiring 9 a y(a, 6, t) — 0. The associated 
points on the surface are shown as a continuous line on 
Fig. ^. The corresponding Eulerian structure is easily 
seen to be an embryonic shock line, parallel to the y-axis, 
with a length O (t 1 ^ 2 ) and a velocity jump also O (t 1 ^ 2 ) , 
except near its end points, where it vanishes. 

It is now easily checked that the separatrix ellipse 
and the zero-Hessian ellipse are tangent at the points 
a = 0, b = ± (-t/(26<)) 1/2 , denoted A and A' on Fig. |. 
These points of vanishing Hessian, which belong to the 
edge of the regular region, are the only kurtoparabolic 
points in the sense of Section [n| Arbitrarily large densi- 
ties are obtained in their neighborhood. Contrary to the 
one-dimensional case, the condition p^ > p, for large 
p, does not put an upper bound, similar to (p4|), on the 
time r elapsed since t*. Actually, in two (and more) di- 
mensions, kurtoparabolic points persist generically for at 
least a finite time, irrespective of the presence of the ad- 



G 



ditional symmetry assumed in the simple example given 
above. 



B. Kurtoparabolic points in two dimensions 

We recall our definition of a kurtoparabolic point A as 
a point (i) where the Hessian of the Lagrangian potential 
<p(a,b,t) vanishes and (ii) which belongs to the bound- 
ary of the regular part of the convex hull f c (ci, b, i). This 
requires two local constraints: that A be parabolic and 
that the surface defined by <p(a, b, t) should not cross the 
tangent plane at A. It also requires a global constraint, 
namely that A should not be situated below a piece of 
the convex hull, which would correspond to A having 
been absorbed by a mature shock before the current time 
t. The latter condition is automatically satisfied at the 
birth of the first singularity and will then persist for at 
least a finite time. The former can be expressed purely in 
terms of the local properties of the Lagrangian potential. 
For this, let us find the normal form associated to a kur- 
toparabolic point at an arbitrary time t (not necessarily 
close to i*). In what follows the time is purely a pa- 
rameter which will not be written. The vanishing of the 
Hessian requires the vanishing of at least one eigenvalue 
of the Hessian matrix of second derivatives. We assume, 
again without loss of generality, that A is the origin, that 
the potential and its gradient vanish at A, that the van- 
ishing eigenvalue of the Hessian matrix corresponds to 
the a-axis and that the other eigenvalue is —fi < 0. We 
now write the local Taylor expansion of the Lagrangian 
potential. Clearly, we need to include terms up to fourth 
order in a, but we shall see that the relevant order in b 
is two (because b — 0(a 2 )). Hence, we can write 



ip(a, b) ~ — A*"o~ aa3 @ a2 b + Co 4 - 



(37) 



We now require that, at least locally (i.e. for small a 
and 6), the surface defined by ip(a,b) should be below 
its tangent plane at the origin. This amounts just to 
{a, b) < and is equivalent to 



0, C < 0, (3 2 < -2K- 



(38) 



(Note that /3 ^ since its vanishing would correspond to 
having a preshock.) The only "sharp" condition is the 
vanishing of the coefficient a of a 3 . Indeed, for a^O the 
surface crosses its tangent plane at the origin. Hence, 
there are two sharp conditions for the existence of a kur- 
toparabolic point: the vanishing of the Hessian and of 
the coefficient of a 3 . Thus, we expect that kurtoparabolic 
points arc found on codimension two manifolds in space- 
time, that is, time-dependent discrete locations which 
persist for at least a finite time. 

Actually, kurtoparabolic points are generally born with 
the first singularity which is itself such a point, albeit one 



with a higher degree of degeneracy (the coefficient (3 also 
vanishes). Generically, kurtoparabolic points disappear 
at large times when only a network of shocks subsists. 
This is seen, for example, in the study of flame-front 
cracks in Ref. ]3?| ]. 

The typical local aspect of a kurtoparabolic point is 
shown in Fig. 0. 



<P (a,b,t) 




FIG. 5. The Lagrangian potential in two dimensions, in 
the neighborhood of a kurtoparabolic point. The continuous 
line is the separatrix between the regular part and the ruled 
surface of the convex hull. The dotted line corresponds to the 
vanishing of the Jacobian of the Lagrangian map. 

The cl i ► cl symmetry, which is here generic, and 
the convexity of sections by planes of constant a allows a 
straightforward construction of the convex hull of (p (a, b). 
The convex hull contains, as in Section IV A, a ruled 
surface made of segments parallel to the a-axis. The 
separatrix is now a parabola of equation 



, 2C a 



(39) 



shown as a continuous line in Fig. [5]. The line of vanishing 
Hessian 



H(a, b) = -2pf3b - 4 (3/< + P 2 ) a 2 = 



(40) 



is also a parabola (shown as a dotted line) which touches 
the former at the kurtoparabolic point. The correspond- 
ing Eulcrian structure is a shock line with an end point 
of the kind shown in Fig. ^. High densities p^ — 
Po/ H(a,b) are obtained in the neighborhood of the kur- 
toparabolic point. Specifically, the set of regular La- 
grangian point such that p( E ' > p is defined by the fol- 
lowing inequalities 



(3b 
2C 



< a 2 < 



po/p + 2p(3b 
~4(3K + /3 2 )' 



(41) 



Using the Lagrangian map x = — d a tp, y = — db^p it may 
be shown that, in the Eulerian space, the mean density 
in a disk of small radius r around the end of the shock 
line p(r) oc r~ 2 / 3 . 
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It follows from ([H]) that, when p is large a and b are 
restricted to being very close to the kurtoparabolic point. 



C. The -7/2 law in two dimensions 

We now assume random and smooth initial conditions 
and proceed along the same general lines as in Sec- 



tion [II B; in particular, spatial periodicity is assumed 



(here and in the next section), although this is not essen- 
tial. The cumulative probability to have > p is now 
expressed in terms of the joint pdf of all the relevant pa- 
rameters at a kurtoparabolic point (its position x„, and 
the three parameters p, (3 and £) in the normal form: 



P>(p;x,t)= f 

Jp< E ) (x,t)>P 



P4(x*, p,, 0, t) cix* dpd(3d(. 



(42) 



(Additional dependences of the probability density on the 
velocity at the kurtoparabolic point and on the orienta- 
tion of the axis corresponding to the vanishing eigenvalue 
are omitted for simplicity, since they are irrelevant for 
the calculation of the density.) An important difference 
with the one-dimensional case, due to the persistence of 
kurtoparabolic points, is the lack of a argument in jxj. 
The probability to have a kurtoparabolic point at x* may 
still depend on x* (unless the initial condition is homo- 
geneous) and will certainly depend on the time t. With 
Gaussian initial conditions it may be shown that is 
non- vanishing for any t > 0. 

Next, we use the very sharp localization of high- 
densities near the Eulcrian points which are associated 
to the kurtoparabolic points and change from Eulerian 
to (local normal form) Lagrangian coordinates, rewriting 
(11) as 

P > ( j o;x,t)~ / H(a,b,t) P4(x, p, /3,£;i) dadb dpdfl <i£, 
Jv 

(43) 

where the domain T> is the set of (a, 6, p, P, £) such that 
(|l|) holds and H (a, b, t), given by ( |34| ) , is the Jacobian of 
the Lagrangian map. In ( |43| ) the integration over a and 
b can be carried out explicitly. This gives P > (p; x, t) oc 
(Po/p) 5 ^ 2 and, hence, 



p(p;x,t) cx — 
\Po 



-7/2 



00, 



(44) 



which is the two-dimensional -7/2 law. (The constant in 
front of the power law, which involves integrals over p, (3 
and is not written.) 

It is of interest to note that we have obtained ex- 
actly the same scaling as in one dimension. The rea- 
son for this is rather interesting. In one dimension, the 



dominant singularities are preshocks which are isolated 
events in space-time, while in two dimensions they are 
kurtoparabolic points which are persistent in time. Nev- 



ertheless, if we compare the integrals (27) and (43) and 
the conditions on the integration domains (28) and (|4l]), 
we find that, in two dimensions the spatial b- variable 
plays exactly the same role as the temporal r-variable in 
one dimension. We can see this also by examining Fig. |^: 
cuts at constant b will have the same shape as the curves 
shown in Fig. [?] for the one-dimensional case, changing 
from a single maximum for b < 0, to a quartically de- 
generate maximum for 6 = and to two symmetrical 
maxima for b > 0. 

We must stress that in two (and more) dimensions, 
the dominant contribution to the tail of the density does 
not come from preshocks as in one dimension but comes 
from the entire life span of kurtoparabolic points which 
are just born at preshocks. Actually, the contribution of 
a small time interval straddling a preshock gives a p" A 
intermediate asymptotic law, going over to a p~'' 7 1" 2 law 
at very large p. 



D. Higher dimensions 

In dimensions d > 2, kurtoparabolic points are now 
space-time manifolds of codimension two. The corre- 
sponding normal form (in suitable coordinates) can be 
written, using a for the coordinate in the direction of the 
vanishing eigenvalue of the Hessian matrix and b = b, t 
(i = 1, . . . , d — 1) for the d — 1 remaining coordinates. 
The analog of (J37|) with a vanishing aa 3 term is now: 



d-l 



^(a,b) ~Ca 4 + E 



A: a h 



where 



d-l 



C < 0, m > 



E-<- 2 c 



(45) 



(46) 



are the analogs of (pq). An easy calculation shows that 
the Jacobian J = (—l) d H(a,h) of the Lagrangian map 
is given by 



/d-l 

J{a, b) = - JJ Pi 



4(3C + E^)« 2 + 2^ / 3 I 6 l 

V i=i ^ i j i=i 



(47) 



Because of the a — > — a symmetry and the convexity in 
all the other variables, the separatrix of the convex hull 
is again obtained by setting d a <p(a, b) = 0, thereby ob- 
taining 



(48) 
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By proceeding as in two dimensions, we can write the 
analog of ( f43| ) in which b, p and j3 must now be reinter- 
preted as (d — l)-dimensional vectors. As to the integra- 
tion domain P, it is now defined by J(a,b) < po/p and 
the negativity of the l.h.s. of (Q), which expresses the 
belonging to the regular domain. 

We now observe that J2i=i Pi^i plays here the same 
role as (3b in the two-dimensional case. There is no b- 
dependence other than this. Hence, assuming that not 
all the /?i's vanish (otherwise we have a preshock) we 
can change variables in the b-space, taking one axis in 
the direction of 53j=i Pik. The integration in all the 
b-directions orthogonal to this direction is trivial and 
gives order unity contributions. The remaining integral 
is just the same as in two dimensions. Hence, in any 
dimension d > 2, the contribution stemming from such 
kurtoparabolic manifolds of codimension two is again a 
p~ 7 / 2 tail in the pdf of the density. 

For d — 3 another type of singularity with vanishing 
Hessian has been identified, denoted by A\ A-$ in Ref. |j4j . 
ft corresponds to a kurtoparabolic point at which the 
tangent (hyper)plane has another point of tangency with 
the graph of the Lagrangian potential. Such points do 
not contribute to the leading order of the pdf for large 
densities. For d > 3, higher-order singularities such as 
the As-points of Ref. |34| can appear generically and we 
do not know how they affect the -7/2 law. 



V. CONCLUDING REMARKS 

In this paper we determined the pdf of the mass den- 
sity for the Burgers equation in the inviscid limit for 
smooth initial conditions which are random but not nec- 
essarily homogeneous. We showed that in one, two and 
three dimensions this pdf has a power-law tail with ex- 
ponent -7/2. In one dimension this tail originates from 
preshocks, that is nascent shocks which take place at dis- 
crete times, whereas in two and three dimensions the tail 
comes from time-persistent boundaries of shocks (associ- 
ated in Lagrangian space to kurtoparabolic singularities). 
In the neighborhood of such points, arbitrary large but 
finite densities are present. Note that, in any case, the 
-7/2 law is due to a phenomenon of shock germination, 
either in space-time or just in space. 

The density pdf will of course be modified if, instead 
of the limit of vanishing viscosity, we assume a finite but 
small viscosity. A similar question has already been con- 
sidered by Gotoh and Kraichnan [^5| for the pdf of large 
(negative) velocity gradients £ in one-dimensional forced 
Burgers turbulence; they found that, at very large values 
of £ there is a power-law range with exponent - I, different 
from the exponent prevailing at those values of £ where 
the inviscid limit is achieved. Similarly, we expect that, 
due to the large shear inside shocks, the -7/2 tail in the 



density pdf would become just an intermediate asymp- 
totic range, beyond which another law should prevail. 




FIG. 6. Position-velocity phase space for the 
(Jeans)-Vlasov-Poisson equation in one dimension. Support 
of the particle distribution at various times for single-speed 
initial conditions. 

As mentioned in the Introduction, the Burgers equa- 
tion with vanishing viscosity is used by cosmologists, un- 
der the name of adhesion model, to approach the prob- 
lem of the formation of large scale structures. In prin- 
ciple, the appropriate mathematical framework should 
involve partial differential equations for density and ve- 
locity, coupled to the Poisson equation for the gravi- 
tational potential. The initial velocity (in an expand- 
ing universe at decoupling) is then determined uniquely 
in terms of initial density fluctuations (see, e.g., Sec- 
tion 2.2.2 of Ref. fl5|l ). When dealing with dustlike colli- 
sionless matter, a hydrodynamical description cannot be 
justified on the usual grounds that local thermodynamic 
equilibrium is quickly achieved. However, as long as par- 
ticles do not cross, a quasi-hydrodynamical description, 
without any viscous diffusion term, is suitable. After 
crossing, the multi-stream situation may be described 
pseudo-microscopically in terms of a distribution func- 
tion /(x, v, i) in the position-velocity phase space satis- 
fying the ci-dimensional (Jeans)-Vlasov-Poisson [p8fl 



a t /(x, v, t) + (v • Vx - Vx • v v ) / = o, 



(49) 



supplemented by the Poisson equation, relating the grav- 
itational potential to the density p(x, t) = J /(x, v, t) dv. 
(We omitted the expansion factor for simplicity.) In 
(|49|), the position and velocity variables are in prin- 
ciple independent, but the relevant solutions for the 
pseudo-microscopic description are of the "single-speed" 
type. By this we understand that the distribution / has 
its support on a d-dimensional submanifold of the 2d- 
dimensional phase space such that, initially, there is a 
single velocity u (x) associated to a given position x. 
Such single-speed solutions may, after particle crossing, 
possess more than one velocity for a given x. The distri- 
bution / does however remain well-defined (see Fig. 
which corresponds to the one-dimensional case). For 
smooth initial data, the support of / remains smooth 
even after particle crossing (see Ref. |39) for the one- 
dimensional case). In one dimension, it is easy to show 
that, near a point x„ where the tangent to the graph of 
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the support is parallel to the w-axis (such as points x± on 
Fig. ||), the density has a singularity oc \x — a;*)" 1 / 2 
When x* is random with a probability density at = 
we infer by an argument similar to that of Section HI B 
that the pdf p(p;x,t) oc p~ 3 for p — > oo. This argument 
carries over to higher dimensions; the a;*-pomts are then 
on the (d — l)-dimcnsional manifold where the Jacobian 
of the Lagrangian map vanishes. 

The same p~ 3 law holds within the (i-dimensional 
Zel'dovich approximation Q . The latter is indeed equiv- 
alent to a modified prescription for determining the grav- 
itational potential; it is obtained from the Poisson equa- 
tion only initially and its gradient (the gravitational 
force) is then taken constant along particle trajectories. 
Clearly, the modified prescription has the only effect that 
it changes the precise position of the support (after the 
first crossing) but not the nature of the ensuing den- 
sity singularities. As noticed in Ref. the power law 
pdf with exponent —3 relative to the velocity gradient 
for the inviscid one-dimensional forced Burgers equation, 
proposed in Ref. (2f| , can be interpreted in terms of sin- 
gularities if one allows multi- valued solutions; this is pre- 
cisely the case in collisionless physical situations. 

The Burgers/adhesion model has been found useful for 
describing large-scale features of collisionless dynamics, 
such as the positions and slow thickening of pancakes. It 
is nevertheless intrinsically a hydrodynamical model and, 
as such, better suited to handle low-pressure collision- 
dominated matter. As far as their singularities are 
concerned, the Zel'dovich approximation and the Burg- 
ers/adhesion model are in different universality classes 
and have different tail behavior for the density pdf. We 
finally mention that with the three-dimensional adhesion 
model, various other quantities may be calculated ana- 
lytically, which are more directly related to what can be 
determined by cosmological observations and/or by N- 
body simulations, such as the density correlation func- 
tion (at small distances) and the pdf of the mass in a 
ball of small radius. 
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